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INTRODUCTION 


A  technique  for  predicting  performance  of  explosive  fragmentation  munitions  presented  in 
this  work  is  based  on  integrating  the  CALE  computer  program  (ref.  1)  with  a  newly  developed 
fragmentation  computer  code  MOTT.  CALE  is  a  two-dimensional  and  three-dimensional  axial 
symmetric  high  rate  finite  difference  computer  program  based  on  Arbitrary  Lagrangian-Eulerian 
formulation  of  the  governing  equations.  The  mathematical  description  of  the  MOTT  code  is 
given  in  a  section  that  follows. 

The  geometries  of  two  example  problems  are  shown  in  figure  1 .  Upon  initiation  of  the  high 
explosive  charges,  rapid  expansion  of  high-pressure  high-velocity  detonation  products  results  in 
high-strain  high-strain-rate  dilation  of  the  hardened  steel  shells,  which  eventually  rupture 
generating  a  “spray”  of  high-velocity  steel  fragments.  In  the  case  of  the  Charge  A,  the  dilation  of 
the  steel  shell  is  accompanied  by  the  implosion  of  the  copper  shaped  charge  liner  that  produces 
a  high-speed  metal  jet  moving  along  the  charge’s  axis  of  symmetry  z.  In  the  case  of  the  Charge 
B,  in  addition  to  the  explosive  charge,  the  hardened  steel  shell  also  encapsulates  a  tracer  and  a 
fuze  that  occupy  significant  volume  of  the  munition  payload.  A  threaded  connection  between 
the  fuse  and  the  fragmenting  portion  of  the  hardened  steel  case  was  also  included  into  the 
CALE  model.  In  order  to  allow  slippage,  the  threaded  joint  was  modeled  with  a  few  rows  of 
computation  cells  employing  ideally-plastic-zero-yield-strength  material  with  the  same 
hydrodynamic  response  parameters  as  steel.  Following  the  expansion  of  the  detonation 
products,  the  fuze  section  of  the  warhead  is  projected  in  the  negative  direction  of  the  z-axis, 
without  contributing  to  the  lethality  of  the  fragment  spray. 

Charge  A  Charge  B 


r  .  Hardened  steel  shell 
Explosive,  A  y 


Threaded  joint  Hardened  steel  shell 
R 


t=0  psec,  V/V0=  1 

Copper  shaped  charge  liner 


t=0  psec,  V/V0= 1 


t=8  psec,  V/V0= 3 


Steel  fragments 


Explosive 


Tracer  cavity 
Fuze 


t  1 3  psec,  V/V0  3  Steel  fragments 

Copper  jet  v 

R 


t=20  psec,  V/V0=14 


t=30  psec,  V/V0=15 


Figure  1 

Results  of  CALE-code  modeling:  initial  configuration  and  CALE’s  predictions  following  the 

explosive  detonation  initiation 
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In  addition  to  the  specification  of  the  problem  geometry  and  initial  and  boundary 
conditions,  equations  of  states  and  constitutive  equations  for  all  materials  have  to  be  specified 
before  the  solution  procedure  can  be  initiated.  The  explosive  was  modeled  using  the  Jones- 
Wilkins-Lee-Baker  equation  of  state  employing  a  set  of  parameters  resulting  from  thermo¬ 
chemical  equilibrium  analyses  of  detonation  products  with  the  JAGUAR  code  and  calibrated  with 
copper  cylinder  test  expansion  data.  The  hydrodynamic  responses  of  the  steel  shell  and  the 
copper  liner  were  modeled  using  a  standard  linear  polynomial  approximation  usually  employed 
for  metals.  The  constitutive  behavior  of  these  metals  was  modeled  using  the  Steinberg-Guinan 
yield-strength  model  and  the  von  Mises  yielding  criterion.  A  standard  set  of  parameters 
available  from  Tipton  (ref.  2)  were  employed  in  the  analyses. 

Since  the  extent  of  dilation  of  the  rapidly  expanding  steel  shell  is  limited  by  its  strength,  at 
some  point  the  shell  ruptures  generating  a  spray  of  steel  fragments  moving  with  trajectories  at 
angles  0  with  z-axis.  Accordingly,  the  principal  topic  of  this  work  is  a  numerical  model  for  the 
analytical  description  of  parameters  of  the  resulting  fragment  spray  as  functions  of  the  “spray” 
angle  0.  In  typical  large-scale  explosive  fragmentation  tests  (arena  tests),  the  tested  munitions 
are  positioned  at  the  origin  of  the  reference  polar  coordinate  system  and  surrounded  with  series 
of  velocity-measuring  screens  and  fragment-catching  witness  panels,  all  at  significant  distances 
from  the  warhead.  Accordingly,  the  fragmentation  characteristics  are  assessed  as  functions  of 
polar  angles  0  'identifying  angular  positions  of  these  measuring  devices.  Assuming  that  the 
fragment  trajectory  angles  0  do  not  change  with  time  (that  is  the  lateral  drift  of  fragments  due  to 
air  resistance  is  small)  and  that  definitions  of  angles  0  and  0  'are  approximately  identical,  the 
developed  model  enables  prediction  of  crucial  characteristics  of  explosive  fragmenting  munitions 
including  the  number  of  fragments,  the  fragment  size  distribution,  and  the  average  fragment 
velocities. 


MOTT  CODE  FRAGMENTATION  MODEL 

For  a  large  part,  the  MOTT  code  fragmentation  model  is  based  on  the  Mott’s  theory  of 
break-up  of  cylindrical  “ring-bombs”  (ref.  3),  in  which  the  average  length  of  the  resulting 
circumferential  fragments  is  a  function  of  the  radius  and  velocity  of  the  ring  at  the  moment  of 
break-up,  and  the  mechanical  properties  of  the  metal.  A  brief  review  of  the  Mott’s  theory  is 
given  here  for  completion.  Following  Mott  and  Linfoot  (ref.  4),  the  “random  variations”  in 
fragment  sizes  are  accounted  through  the  following  fragment  distribution  relationship 

b/ v/2 

x  ,,  (1) 

In  equation  1 ,  N(m)  represents  total  number  of  fragments  of  mass  greater  thanm ;  p  is 
defined  as  one  half  of  the  average  fragment  mass,  N0=M  ///and  M  is  the  total  mass  of  the 
fragments. 

In  attempting  to  evaluate  the  distribution  of  fragment  sizes  occurring  in  the  dynamic 
fragmentation  of  expanding  metal  rings,  Mott  (ref.  3)  introduced  an  idealized  model  in  which  the 
average  circumferential  fragment  lengths  are  not  random,  but  determined  by  the  interaction  of 
stress  release  waves  originating  from  instantaneous  fractures  in  the  body.  A  schematic  of  the 
Mott’s  model  is  shown  in  figure  2  (a).  Assuming  that  a  fracture  in  the  ring  is  supposed  to  have 
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occurred  first  at  A 1  and  that  stress  release  waves  have  traveled  to  points  B,  and  Bh  further 
fractures  can  no  longer  take  place  in  regions  A1B1  and  A^BV  On  the  other  hand,  in  the  regions 
and  B ,£>  the  plastic  strain  is  increasing,  which  increases  the  probability  of  fractures  at  any 
point  in  these  regions,  especially  at  points  Bu  Bly  £,,  and  Thus,  according  to  Mott’s  theory, 
the  average  size  of  fragments  is  determined  by  the  rate  at  which  stress  relieved  regions  A ^ 
and  A-\B\  spread  through  the  plastically  expanding  ring. 


Stress  Release  Wave 


Figure  2 

Fragmentation  model 

At  the  moment  of  fracture,  let  r  be  the  radius  of  the  ring  and  V  be  the  velocity  with  which 
the  shell  is  moving  outwards.  Then,  according  to  Mott  (ref.  3),  the  average  circumferential 
length  of  the  resulting  fragments  is 


*n  = 


'2P^m 
PY 


(2) 


In  equation  2,  p  and  PF  denote  the  density  and  the  strength,  respectively;  and  y  is  a  semi- 
empirical  statistical  constant  determining  the  dynamic  fracture  properties  of  the  material. 
Given  that  the  shape  and  the  average  fragment  lengths  are  known,  the  idealized  average 
fragment  mass  can  be  calculated.  For  example,  assuming  approximate  cubic-shaped 
fragments,  p  takes  the  following  form 


A  schematic  for  the  newly  developed  technique  implemented  in  the  MOTT  code  is  shown 
in  figure  2  (b).  For  computational  purposes,  the  shell  is  discretized  into  a  finite  number  of  short 
“ring"  segments,  N  .  For  each  discrete  ring  element,  j  uniform  field  variables  are  assumed. 
Accordingly,  the  masses,  the  velocities,  and  radii  of  ring  segments  j  are  defined  by  the  mass 
averages  of  the  respective  parameters 
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(4) 


mi = Zm- 


Z  R>m> 

_h _ L_ 

mj  sin0; 


(5) 

(6) 


0.  <0,.  <0,+  — 
1  2  N  '  1  2  N 


(7) 


In  equations  4  through  6  m, ,  v( ,  and  /?,  denote  the  mass,  the  velocity,  and  the  radial  coordinate 
of  the  /-/^computational  cell  from  the  CALE-code  generated  data.  L;  denotes  a  number  of 
computational  cells  contained  in  the  j-th  ring  segment.  0;  denotes  the  0-angle  that 
corresponds  to  the  j-th  ring  segment  given  by 


0,= 


(8) 


For  each  computational  cell  i ,  the  velocity  v,  and  the  0  angle,  0, ,  are  calculated,  respectively, 
by 


v,=^l 


+  v 


and 


0,  =  arctan  — 


(9) 

(10) 


In  equations  9  and  10,  vzi  and  vRi  denote  the  axial  and  the  radial  velocity  components  from  the 
CALE-code  generated  data. 

Given  that  the  velocities  and  the  radii  of  ring  segments  j  are  determined  through 
equations  5  and  6,  the  resulting  fragment  size  distributions  in  each  segment  j  can  be  calculated 
through  the  following  relationships 


Nj(m)  =  N0Je 


(11) 


hi  = 


I1 

f pF ^ 

3/2 

M 

ip 

7  ) 

lvi) 

(12) 
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(13) 


N  .  = 

°J  Mj 


As  the  detonation  wave  travels  along  the  shell  length  and  the  expanding  detonation 
products  rupture  the  shell,  in  the  case  of  the  idealized  “long-pipe-bombs,”  the  break-up  radii  r; 

and  the  break-up  velocities  V}  of  the  individual  segments  j  are  approximately  the  same, 
regardless  of  the  axial  positions  of  the  segments.  Accordingly,  taking  p=  gj ,  the  number  of 
fragment  distribution  relationship  is  given  by  the  original  Mott’s  equation  1 . 

However,  in  the  case  of  conventional  explosive  fragmentation  munitions  with  shell 
geometries  far  from  that  of  the  idealized  “long-pipe-bombs,”  the  break-up  radii  rj  and  the  break¬ 
up  velocities  V.  vary  along  the  shell  length,  so  that  the  resulting  variance  in  the  average 
fragment  half-weights  of  the  individual  segments  j  may  be  rather  significant.  The  existence 

of  significant  differences  in  the  average  fragment  sizes  between  the  cylindrical  and  the  curved 
portions  of  the  shell  had  been  experimentally  confirmed  in  this  work  through  flash  radiography 
and  high-speed  photography.  Accordingly,  the  following  two  fragment  distribution  relationships 
are  introduced  herein.  The  “shell-averaged”  fragment  distribution  is  defined  as 


N(m)=N/^ 

(14) 

where  Nu  and  are  defined  as 

j 

(15) 

=  IN 

1  0 

(16) 

The  “ring-segment-averaged”  fragment  distribution  is  defined  as 

j 


(17) 


MOTT  CODE  VALIDATION:  CHARGE  A  ANALYSES 

t 

The  validation  of  the  MOTT  code  fragmentation  model  was  accomplished  using  the 
existing  Charge  A  arena  test  data.  The  fundamental  assumption  of  all  fragmentation  analyses 
presented  in  this  work  was  that  the  fragmentation  occurs  instantly  throughout  the  entire  body  of 
the  shell.  Following  Mott’s  critical  fracture  strain  concept  (ref.  3)  and  assuming  that  for  given 
shell  geometry  and  materials,  the  shell  fragmentation  time  is  a  function  of  the  cumulative 
dilatational  plastic  strain  in  the  shell,  the  shell  fragmentation  time  can  be  conveniently 
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expressed  in  terms  of  the  global  shell  dilatational  properties.  Given  that  in  a  typical  fragmenta¬ 
tion  munition  device  the  explosive  is  tightly  confined  inside  the  shell,  the  cumulative  strains  of 
the  expanding  explosive  and  the  surrounding  shell  are  nearly  proportional.  Accordingly,  the 
critical  fracture  strain  at  the  moment  of  the  shell  break-up  may  be  conveniently  measured  in 
terms  of  the  high  explosive  detonation  products  volume  expansions,  V/V0. 

Figure  3  shows  the  effect  of  the  shell  fragmentation  time  on  the  fragment  spray  velocity 
distribution  function.  The  seemingly  significant  disagreement  betweeh  the  experimental 
velocities  and  the  analyses  for  0<1 5°  is  due  to  deliberate  omission  of  the  shaped  charge  jet  data 
from  the  MOTT  fragmentation  analyses,  mainly  because  of  the  minimal  contribution  to  the 
overall  fragment-spray  lethality.  Accordingly,  the  copper  shaped  charge  jet  was  neglected  in  all 
fragmentation  analyses,  although  included  in  the  CALE  model  in  order  to  maintain  proper 
explosive  confinement  parameters.  As  shown  in  figure  3,  varying  the  shell  fragmentation  time 
from  approximately  8  ps  (at  which  the  detonation  products  had  expanded  approximately  three 
times  its  original  volume,  WV0~3  to  approximately  20  ps  (V7V0~14),  the  resulting  changes  in 

the  fragment  spray  angles  0  were  rather  small,  while  the  fragment  spray  velocities  were  affected 
rather  significantly.  As  shown  in  figure  3,  delaying  the  moment  of  the  shell  break-up  predicts 
considerable  increases  in  the  fragment  spray  velocities,  apparently  due  to  the  prolonged 
“pressurized”  interaction  with  expanding  detonation  products  that  increased  the  total  momentum 
transferred  to  the  shell. 


Figure  3 

Fragment  velocity  distribution  versus  spray  angle  0  change  A 

Figures  4  through  6  show  plots  of  the  number  of  fragments  in  the  fragment  spray  as 
functions  of  the  fragment  size  ml /j  ,  the  choice  of  the  number  of  fragment  distribution  function 
model,  equations  14  and  17,  the  spray  angle  0,  the  shell  fragmentation  time,  and  the  dynamic 
fracture  parameter  y.  As  shown  in  the  figures,  increases  in  the  parameter  y  had  resulted  in 
increases  of  the  number  of  fragments  N  ,  both  for  the  N-ml/u  and  the  N-Q  relationships. 
These  results  are  in  agreement  with  the  Mott’s  theory  (ref.  3),  according  to  which  the  parameter 
y  defines  the  probability  of  fracture  in  the  plastically  expanding  shell  determining  the  number  of 
breaks  in  the  circumferential  direction. 
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Figure  4  shows  a  plot  of  a  series  of  curves  given  by  equation  14,  N  {m)= f(X4)(m,y) ,  all 

analyses  repeated  for  two  parameters  considered:  the  shell  fragmentation  time  assumed  and 
the  dynamic  fracture  constant,  y.  For  example,  taking  the  8  ps  (V/Vo-3.0)  fragmentation  time 

with  y=12  and  the  20  ps  (V /V0~14)  fragmentation  time  with  y=30  resulted  in  nearly  identical 

fragment  distribution  curves,  both  in  good  agreement  with  the  data.  The  accepted  shell 
fragmentation  time  was  determined  from  the  high-speed  photographic  data  of  Pearson  (ref.  5). 
Following  Pearson  (ref.  5),  the  fragmentation  of  shells  with  the  idealized  cylindrical  geometries 
occurs  approximately  at  3  volume  expansions,  the  instant  of  fragmentation  defined  as  the  time 
at  which  the  detonation  products  first  appear  emanating  from  the  fractures  in  the  shell. 
Accordingly,  the  accepted  shell  fragmentation  time  was  approximately  8  ps  ( V/V0~3.0 ). 
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Cumulative  number  of  fragments  in  the  fragment  spray  versus  the  fragment  size  ml/j  charge  A 

Figure  5  shows  a  plot  of  the  cumulative  number  of  fragments  from  the  arena  test  data  and 
the  CALE-MOTT  analytic  predictions  employing  N(m)=f(u)(m,y)  and  N(m)= f(xl){m,y) 

relationships;  i.e.,  the  overall  “shell-averaged”  fragment  size  distribution,  equation  14,  and  the 
"ring-segment-averaged”  fragment  size  distribution,  equation  17,  respectively.  As  shown  in  the 
plot,  employing  equation  14  resulted  in  consistently  higher  values  of  the  cumulative  numbers  of 
fragments  than  that  of  equation  1 7,  apparently  because  of  the  nature  of  the  definition  of  the 
overall  “shell-average”  fragment  mass  p0 ,  equation  16.  As  shown  in  the  figure,  two  equally 

reasonable  fits  were  obtained  for  both  relationships  considered,  resulting  in  y=12  for  equation  14 
and  in  y=14  for  equation  17  with  the  standard  deviations  of  approximately  cry=n  (l.62)=2%  and 

cry-l4  (l  .62)=2% ,  respectively.  Both  curves  were  fitted  at  a  single  point  ml fj=  1 .62 ,  which 

corresponds  to  the  total  number  of  fragments  with  mass  greater  than  3  grains;  the  aim  was  to 
replicate  the  overall  lethality  of  the  fragmenting  spray,  rather  than  focusing  on  the  entire  range  of 
the  ml ju  values  considered.  Accordingly,  the  accepted  shell  fragmentation  had  been 

approximately  8  ps  ( V/V0~3.0)  and  two  values  y=12  and  y=14  were  selected  for  all  further 
analyses. 
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Figure  5 

Cumulative  number  of  fragments  versus  ml [i  for  varying  y  CALE-MOTT  analyses  using 

equations  14  and  17  -  change  A 

Figure  6  shows  a  plot  of  the  number  of  fragments  with  mass  greater  than  3  grains  versus 
the  spray  angle  0,  which  is  the  principal  lethality  parameter  of  the  fragment  spray  of  the 
munition.  The  disagreement  between  the  analyses  and  small  spikes  at  45  and  60  deg  is  due  to 
fragments  from  the  shaped  charge  liner-retaining  ring,  which  had  not  been  included  in  the  CALE 
model,  mainly  because  of  the  minimal  effect  on  the  overall  fragment  lethality.  The  disagreement 
between  the  analyses  and  the  spike  at  155  deg  is  due  to  fragments  from  a  rotating  band  that 
had  not  been  included  in  the  CALE  model.  As  shown  in  the  figure,  even  using  a  relatively  crude 
assumption  of  the  shell  fragmentation  time,  the  overall  agreement  between  the  analyses  and  the 
experimental  data  is  very  good. 
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CHARGE  B  MODELING  AND  EXPERIMENTATION 

After  having  established  the  crucial  parameters  of  the  model,  Charge  B shown  in  figure  1 
was  designed  and  optimized.  Upon  fabrication,  the  performance  of  the  new  charge  was  tested 
in  a  series  of  experiments  including  flash  radiography,  high-speed  photography,  and  sawdust 
fragment  recovery. 

The  flash  radiography  tests  were  performed  using  two  150  kV  x-ray  heads  located 
approximately  74  in.  in  front  of  the  round.  Shortly  after  initiating  the  round,  each  of  the  two  x- 
rays  heads  were  flashed  at  the  separate  prescribed  times  and  intervals  several  microseconds 
apart.  Two  flash  radiography  tests  were  conducted.  Each  test  resulted  in  two  dynamic  images 
of  the  expanding  fragmented  steel  shell,  both  images  superimposed  on  the  film. 

The  high-speed  photography  tests  were  performed  employing  a  Cordin  Framing  Camera, 
Model  No.  121 ,  capable  of  recording  up  to  26  high-speed  exposure  frames  with  time  intervals 
between  individual  frames  of  less  than  1  ps  apart.  In  the  experiments,  the  round  was  placed  on 
a  test  stand  in  front  of  a  fiducial  grid,  surrounded  with  four  Argon  gas  light  bombs,  all  enclosed  in 
a  white  paper  tent.  A  total  of  two  high-speed  photography  tests  were  conducted,  each  test 
resulting  in  over  20  dynamic  images  of  the  expanding  and  fragmenting  shell,  approximately  1  us 
apart. 


Figure  7  shows  a  comparison  between  the  CALE  code  predictions  and  the  images  of  the 
expanding  and  partially  fractured  shell  obtained  from  the  flash  radiography  and  high-speed 
photography  experiments.  The  figure  shows  that  the  model  resulted  in  an  accurate  prediction  of 
the  shape  of  the  expanding  hardened  steel  shell,  including  the  early  break  out  of  the  detonation 
products  through  the  joint  between  the  fuze  and  the  main  charge. 
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Early  break-up 


Detonation  products  break  out  through 
cracks 


Figure  7 

CALE  code  modeling  and  experimentation  -  charge  B 
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Referring  to  the  flash  radiographic  image  of  the  partially  fragmented  shell  given  in  figure  7, 
it  is  important  to  note  a  remarkable  difference  between  shapes  and  sizes  of  fragments  ejected 
from  the  cylindrical  and  the  curved  portions  of  the  charge.  As  seen  from  both  the  27  ps  and 
from  the  45  ps  radiographic  images,  the  majority  of  cracks  in  the  cylindrical  portion  of  the  charge 
are  in  the  axial  direction,  resulting  in  a  fragment  spray  with  relatively  large  axially  oriented 
splinter-like  fragments.  On  the  contrary,  in  the  curved  nose  portion  of  the  shell,  the  orientation 
of  cracks  is  random  and  the  distances  between  fissures  are  shorter,  resulting  in  a  fragment 
spray  of  predominantly  small  compactly  shaped  fragments. 

Figure  8  shows  a  series  of  high-speed  photographic  images  of  the  condition  of  the  surface 
as  the  hardened  steel  shell  expands  and  the  detonation  products  emerge  through  the  cracks. 

As  shown  in  the  figure,  visible  fractures  start  to  appear  on  the  surface  of  the  shell  at 
approximately  9.4  ps  after  the  detonation,  which  according  to  the  CALE  code  analyses 
corresponds  to  approximately  V/V0 =1 .8.  Examination  of  the  entire  series  of  images  taken 

approximately  1  ps  apart  shows  that  as  the  shell  expands,  fractures  develop  first  in  the 
cylindrical  portion  of  the  shell.  As  the  shell  continues  to  expand,  the  developed  fractures 
propagate  mostly  in  the  axial  direction,  occasionally  linking-up  by  new  cross-cracks  in  the 
circumferential  direction,  ultimately  resulting  in  formation  of  large  splinter-like  fragments  shown 
in  the  radiographic  images  of  figure  7. 
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Development  of  surface  fractures  in  the  expending  shell  -  charge  B 


As  the  shell  continues  to  expand,  fractures  gradually  advance  towards  the  curved  nosed 
portion  of  the  shell,  and  at  approximately  19.8  ps,  or  at  approximately  6.2  volume  expansions, 
the  entire  shell  is  fully  fragmented,  the  fragmentation  being  defined  as  the  instant  at  which  the 
detonation  products  first  appear  as  they  emanate  from  the  fractures  in  the  shell.  Following  the 
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CALE-MOTT  model  assumption  that  the  critical  fracture  strain  at  the  moment  of  the  shell  break¬ 
up  is  expressed  in  terms  of  the  high  explosive  detonation  products  volume  expansions,  the 
“average”  volume  expansion  at  the  time  of  the  shell  break-up  is  then  approximately  one-half  of 
the  value  of  volume  expansions  of  the  fully  fragmented  shell,  hence  V/V0  =3.1 .  It  is  interesting 

to  note  that  the  value  of  V/V0=3.1  is  in  excellent  agreement  with  that  assumed  initially  based  on 
the  high-speed  photographic  data  of  Pearson  (ref.  5). 

Figure  9  shows  plots  of  CALE-MOTT  analytic  predictions  of  the  fragment  velocity 
distribution  function  and  the  experimental  data.  The  experimental  data  considered  here  was 
reduced  from  the  radiographic  images  of  the  expanding  and  partially  fractured  shell.  As  shown 
in  the  figure,  two  analytical  fragment  velocity  distribution  functions  were  considered:  (i)  assuming 
that  the  entire  shell  fragments  instantly  at  approximately  at  13  ps  ( V/V0=3),  and  (ii)  assuming 

that  the  entire  shell  fragments  instantly  at  approximately  at  30  ps  ( V /V0  =15).  Given  that  the 

explosive  was  modeled  using  a  semi-empirical  set  of  parameters  calibrated  with  the 
experimental  copper  cylinder  expansion  data,  the  CALE  predictions  of  the  expanding  (but  not 
fractured)  shell  velocities  should  be  reasonably  accurate.  As  shown  in  the  figure,  the  agreement 
between  the  data  and  the  V /V0  =3  curve  is  significantly  better  than  with  the  V/V0  =15  curve, 

suggesting  that  the  shell  fractured  at  approximately  3  volume  expansions.  In  addition,  given  that 
(according  to  the  high-speed  photography)  the  entire  shell  had  fully  fractured  at  approximately 
6.2  volume  expansions,  the  assumption  of  the  V/V0=3  instantaneous  fracture  event  is  quite 

reasonable:  once  the  detonation  products  start  to  escape  through  the  cracks,  any  further  gains 
in  the  velocities  are  relatively  small  and  may  be  neglected. 


Figure  9 

Fragment  velocity  distribution  versus  spray  angle  0  -  charge  B 
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The  fragment  recovery  tests  considered  in  this  work  were  performed  employing 
disposable  plastic  containers  measuring  approximately  52  in.  in  diameter  and  59  in.  in  height, 
filled  with  approximately  1 000  lbs  of  sawdust.  After  inserting  the  tested  round  in  an  inflatable 
rubber  balloon  and  filling  the  balloon  with  air,  the  balloon  was  positioned  in  the  sawdust  in  the 
middle  of  the  plastic  container.  After  detonating  the  round,  the  resulting  fragments  were 
recovered  employing  a  combination  of  the  magnetic  (for  separating  fragments  from  the  sawdust) 
and  the  vacuum  (for  separating  sawdust  from  the  fragments)  recovery  techniques.  The  mass 
distribution  of  the  collected  fragments  was  analyzed  employing  an  electronic  high  precision 
balance  gauge  interfaced  with  a  computer  system  capable  of  automatic  counting  of  fragments 
as  they  were  placed  on  the  scale  and  weighed.  Weighing  of  the  fragments  was  performed 
employing  the  Ohaus  Voyager  Balance  model  VI 4 130  gauge  with  the  maximum  capacity  of  410 
g  and  the  precision  of  0.001  g.  A  total  of  two  sawdust  recovery  tests  were  performed,  each  of 
the  tests  resulting  in  successful  recovery  of  approximately  over  99.8%  of  the  mass  of  the  steel 
shell. 


Figure  1 0  shows  plots  of  CALE-MOTT  analytic  predictions  of  the  cumulative  number  of 
fragments  compared  with  the  data  from  the  fragment  recovery  tests.  As  shown  in  the  figure,  two 
analytic  relationships  had  been  considered:  (i)  the  “shell-averaged”  fragment  size  distribution  (eq 
14)  and  (ii)  the  “ring-segment-averaged”  fragment  size  distribution  (eq  17).  As  shown  in  the 
figure,  the  analytic  prediction  of  equation  14  significantly  disagrees  with  the  experimental  data, 
regardless  of  the  value  y  considered.  The  disagreement  between  the  equation  14  predictions 
and  the  data  is  mostly  because  of  the  significant  variance  in  fragment  weights  //,  along  the 

shell,  ultimately  resulting  in  over-predicting  the  “shell-averaged”  fragment  weight  //„  (eq  16). 

On  the  contrast,  the  agreement  between  the  data  and  the  y=14  “ring-segment-averaged” 
fragment  size  distribution  given  by  the  equation  17  is  quite  good:  <r^l4  (0.1944)=-7.3%.  Given 

relative  simplicity  of  the  model,  the  overall  agreement  between  the  analyses  and  the  data  is 
excellent. 


Figure  10 

Cumulative  number  of  fragments  versus  the  fragment  mass  m  -  charge  B 
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CONCLUSIONS 


A  technique  for  predicting  performance  of  explosive  fragmentation  munitions  presented  in 
this  work  is  based  on  integrating  three-dimensional  axisymmetric  hydrocode  analyses  with 
analyses  from  a  newly  developed  fragmentation  computer  code  MOTT.  The  validation  of  the 
MOTT  code  fragmentation  model  was  accomplished  using  the  existing  munition  arena  test  data. 
After  having  established  the  crucial  parameters  of  the  model,  a  new  explosive  fragmentation 
munition  was  designed  and  optimized.  Upon  fabrication  of  the  developed  munition,  the 
performance  of  the  new  charge  was  tested  in  a  series  of  small-scale  experiments  including  the 
flash  radiography,  the  high-speed  photography,  and  the  sawdust  fragment  recovery. 

Considering  relative  simplicity  of  the  model,  the  accuracy  of  the  MOTT  code  predictions  is  rather 
remarkable. 
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